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@ Method and device for processing a signal. 

(57) A method and device for processing a signal 
are described, wherein an estimate of a mul- 
tipath-induced contribution to a demodulated 
navigation signal is calculated and subtracted 
from said demodulated navigation signal to 
obtain an estimated line of sight contribution to 
said demodulated navigation signal, and a pro- 
pagation time t 0 is calculated from the thus 
calculated line of sight contribution to said 
demodulated navigation signal, such that a very 
accurate propagation time t 0 of said navi- 
gation signal can be calculated. 



c 

z3™ - 






-I £E 

22 

C < 0 
— tc ° 


o □ « 

0)2 

w 
o 


' CLC 


in 




« £ w 

ii Z 
Ul Ul 






a 

10 


CO 



DC 

t2s 

goo 
o 55 o 
c 
a. 



< _ 

CCu> 
0«» 



5« 

EC m 



^5 



<5 



g 

Si 



1:" 



UJ 

a. 



Jouve, 18, rue Saint-Denis, 75001 PARIS 



EP 0 654 677 A1 

BACKGROUND OF THE INVENTION 
1. Field of the invention. 

5 The present invention relates to a method and device for processing a signal, for use in communication 

systems or satellite navigation systems, for instance. More particularly, the present invention relates to a meth- 
od and device for processing a signal wherein any errors caused by multipath are eliminated or at least reduced. 
In one particular embodiment, the present invention relates to a navigation system. 

10 2. Description of the Prior Art. 

In principle, communication from a transmitter to a receiver involves the transmission, by the transmitter, 
of a carrier wave signal which is coded by an information signal, such as by AM or FM, and the reception, by 
the receiver, of the coded signal. In the receiver, the received signal is decoded for obtaining the information 

15 therein. The propagation time of the signal, i.e. the time which passes between the moment of transmission 
and the moment of reception, depends upon the distance between the transmitter and the receiver, and upon 
the propagation speed of the carrier wave signal, as is well known. 

In practice, the carrier wave signal reaches the receiver not only directly, i.e. via the line of sight, but also 
via one or more reflections from objects such as buildings, mountains, clouds, atmospheric layers, etc. Since 

20 the propagation path via reflection is always longer than the line of sight the propagation time of the signal 
via reflection is always longer than the propagation time of the signal via the line of sight, i.e. a reflected signal 
always has a certain delay with respect to the direct signal. The receiver can not distinguish between the direct 
signal and a reflected signal, and the combination of the direct signal as received by the receiver and the re- 
flected signals as received by the receiver is processed as if it were an undistorted signal; however, said com- 

25 bination obviously contains an error with respect to the signal as transmitted. This problem is commonly known 
as "multipath". Multipath can also be caused by diffraction. 

Communication of the above-described kind is commonly used for many applications, and in most appli- 
cations it is a general desire to reduce the errors caused by multipath. However, one application where reducing 
such errors is extremely important is navigation, such as satellite navigation, and the present invention will be 

30 described with reference to a navigation system, although it is explicitly noted that the invention is not restricted 
to such use. 

In the context of the present application, navigation is based on the principle of determining the distance 
between an object and a reference position by means of a communication signal. A reference station, such as 
a satellite, emits a signal, which is received by a receiver associated with said object. By measuring the prop- 

35 agation time of the signal, the distance of the object to the reference station can be calculated. Similarly, the 
distances of the object to a second and possibly further reference stations are calculated. In the case of a sta- 
tionary reference station, the position of the reference station is constant and may be incorporated in the re- 
ceiver as a constant value. In the case of a non-stationary reference station such as a satellite, the signal may 
contain information regarding the position of the reference station. Since, therefore, the positions of the refer- 

40 ence stations are "known" by the receiver, it is thus possible to calculate the exact position of the object. 

Examples of such navigation systems of the above-described type are known under the names GPS and 
LORAN-C. These systems are used, for instance, by ships at sea, and by aeroplanes in their approaching an 
airport. For a more detailed description of such systems, reference is made to the following publications: 
J.J. Spilker: "Digital Communications by Satellite", Prentice Hall, New Jersey, 1977; 

45 N. Ackroyd and R. Lorimen "Global Navigation: A GPS User's Guide", Lloyd's of London Press, London, 1990; 
B. Forsell: "Radio Navigation Systems", Prentice Hall, 1991, ISBN 0-13-751058-6. 

It will be evident that in the above-mentioned applications, accuracy is of utmost importance, and any error 
due to multipath should be eliminated or at least minimized. 

50 SUMMARY OF THE INVENTION 

It is, therefore, a general objective of the present invention to eliminate multipath-induced errors. 

This and other objects are achieved in accordance with the present invention by producing a method and 
a device for processing a received signal, said method and device having the capability of detecting the pres- 
55 ence of multipath and the capability of determining the contribution of the line of sight signal and the contri- 
bution of the multipath signals. The inventive method and device then removes the contribution of the multipath 
signals from the received signal, thus leaving an uncorrupted line of sight signal. As such, using the uncorrupted 
line of sight signal, an accurate propagation time from a signal source to the receiver can be determined. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be explained in more detail by discussing a preferred embodiment with reference to the 
drawing, wherein: 

fig.1 schematically shows a situation for employing the present invention; 
f ig.2 schematically shows the shape of an undistorted navigation signal; 
fig.3 illustrates a correlation operation; 

f ig.4 schematically shows the idealized shape of a correlation function, and illustrates the conventional 

method of processing the navigation signal; 
fig.5 schematically shows a combination of three correlation function peaks of a multipath signal, and 

illustrates the drawbacks of the conventional method; 
f ig.6 illustrates the steps taken by a preferred method of the invention; 

f ig.7A schematically shows the idealized shape of a correlation function, and illustrates the inventive meth- 
od of processing the navigation signal; 
f ig.7B schematically shows the shape of an auxiliary function; 
f ig.8 illustrates the effect of the invention; 

f ig.9 illustrates an example of a device according to the invention; and 
fig. 10 illustrates the shape of a correlation function in practice. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Fig.1 schematically shows a ship 1 at sea, and a satellite 2 emitting a signal 4. The satellite equipment for 
generating and emitting said signal is generally indicated at 3. The construction of the satellite equipment 3 
does not form part of the invention, and knowledge of this construction is not necessary for understanding the 
present invention. Therefore, the satellite 2 and its satellite equipment 3 will not be discussed further. Suffice 
it to say that the satellite 2 may be a satellite as being used at present, such as the GPS satellites. 

For the purpose of the present discussion it will be neglected that the signal 4 contains data carrying an 
information message, and only a navigation signal 5 modulated on the carrier wave of the signal 4 is considered. 
Fig.2 schematically shows the shape of the navigation signal 5 as emitted by the equipment 3. In fig.2, the 
ordinate is the signal intensity I in arbitrary units, and the abscissa is the time t. The navigation signal 5 may 
have the form of a single pulse, such as is employed in the LORAN-C system (curve A), or the navigation signal 
5 may have the form of a sequence of pulses (curve B), such as is employed in the GPS system. Fig.2 illustrates 
that the navigation signal 5 is emitted periodically, the period P being 1 ms in the above-mentioned systems. 
It is noted that the beginning of each signal period is emitted at moments to in time which are very accurately 
specified under control of a system clock. 

Curve A of fig.3 illustrates the GPS navigation signal 5 at a larger time scale, wherein only a first part of 
the signal is shown. The GPS navigation signal 5 constitutes a pseudo random noise code containing a pre- 
determined number of code bits. In order to distinguish the code bits from data bits, the code bits are usually 
referred to as "chips". Each chip can have two possible values, namely +1 or-1. In the GPS system, said pre- 
determined number is 1023. In the following, the individual chips will be referred to as C(1), C(2), ... C(1023). 
Each satellite 2 has its own characteristic code as defined by the set of values of C(1), C(2), ... C(1023). 

The receiver 1 0 aboard the ship 1 comprises at least one code generator, which can be set to generate a 
reference signal corresponding to the characteristic code of a specific satellite 2. In general, the receiver 10 
is equipped with a memory where information on the characteristic codes of a plurality of satellites is stored, 
so that an operator needs only to indicate which satellite he is interested in and the code generator of the re- 
ceiver 1 0 will automatically generate a reference signal corresponding to the characteristic code of that specific 
satellite of interest (tuning the code generator). 

Figures 3 and 4 illustrate the principle of detecting the propagation time of the signal 4, which will be re- 
ferred to as delay time t 0 . As mentioned above, curve A of fig.3 shows the navigation signal 5 as emitted by 
the satellite 2, i.e. the beginning of C(1) coincides with the beginning to of a time period. Therefore, curve A of 
fig.3 can also be considered a representation of the reference signal as generated by the code generator of 
the receiver 10. Curve B shows the navigation signal 5' as received at an input 11 of the receiver 10 (in the 
absence of noise and multipath), i.e. the beginning of C(1) being delayed by delay time t 0 with respect to t©. It 
is noted that in fig.3 the carrier wave and data are not shown. 

In order to determine the time delay t 0 , a control device (not shown) associated with the receiver 10 pro- 
vides a local reference signal 6, which is identical to the reference signal 5 but shifted over a certain shift time 
x. This local reference signal 6 is shown in curve C of fig.3. 

In order to determine whether the (known) shift time x equals the (unknown) time delay t 0 , said control 
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device compares said local reference signal 6 with the navigation signal 5' as received at an input 11 of the 
receiver 10. More specifically, said control device performs a multiplication operation on the local reference 
signal 6 and the received navigation signal 5', i.e. the local reference signal 6 is multiplied point by point with 
the received navigation signal 5'. Subsequently, the resulting signal, which is shown in curve D of fig. 3 and 
5 which will be referred to as multiplied signal 7, is averaged, i.e. the multiplied signal 7 is integrated over a suf- 
ficiently large time to provide a mean value M(t) of the multiplied signal 7. 

The multiplied signal 7 can only take the values +1 or -1, depending on whether the values of the signals 

6 and 5' are identical to each other or not, respectively. Therefore, the mean value M(x) of the multiplied signal 

7 will be a value between +1 and -1 , too. Since the shape of the multiplied signal 7 depends on the shift time 
10 x, the exact value of M(x) depends on the chosen value of the shift time t. Since the navigation signal 5 is a 

pseudo random code, said mean value M(x) will be practically zero for most values of the shift time x, and will 
differ significantly from zero only if the shift time x approaches the delay timex 0 . In an ideal case, if the shift 
time x equals the delay time x 0 , said mean value M(x) will be equal to 1 (one). Curve E of fig. 3 illustrates this 
characteristic of said mean value M(x) as a function of the shift time x, this function being referred to as cor- 
15 relation function 8. 

The exact shape of the correlation function 8 is known in advance, and depends on the specific code of 
the satellite in question. In curve E of fig.3, this shape is shown as being an ideal triangle. The width W of the 
triangle is a known constant value, for practical purposes equal to twice the duration of a chip; therefore, the 
correlation function 8 is fully determined by the coordinates (x 0 , M(x 0 )) of its maximum. Herein x 0 is the relevant 
20 parameter determining the position of the correlation function 8 with respect to the beginning of the time per- 
iods. 

It is noted that in practice the correlation function 8 is not an ideal triangle as shown in curve E of fig.3 
but is more or less curved, as is shown in fig. 10, depending on the characteristics of the receiver. 

Further it is noted that in systems such as LORAN-C, wherein in each time period only one navigation signal 
25 pulse having a predetermined shape is emitted instead of a pseudo random code, the navigation signal pulse 
itself can be substituted for the correlation function 8 as discussed above. 

In the following, the conventional method employed in conventional receivers for determining the propa- 
gation time x 0 is explained with reference to fig .4, which shows the correlation function 8 on a larger scale. 
Firstly, the local reference signal 6 is generated for a first value t 1 of the shift time x, and the corresponding 
30 mean value M(t,) of the multiplied signal 7 is calculated. Secondly, the local reference signal 6 is generated 
for a second value ^+At of the shift time, and the corresponding mean value M^+At) of the multiplied signal 
7 is calculated. Herein, the sample spacing, indicated as At, is a constant value. 

From the respective values M(t,) and M(ti+At) an error signal e, = M(t,) - M(tt+At) is calculated. This error 
signal is used to control the shift time of the local reference signal 6, such that when the mean value M(x) is 
35 subsequently calculated for two values t 2 and t 2 +At, the error signal e 2 = M(y - M(t 2 +At) is less than e,. This 
procedure, which is known as Delay Locked Loop (DLL), is executed continuously so that in practice at any 
particular moment e « 0, and, based on the assumption that the peak of the correlation function 8 is symmet- 
rical, the propagation time x 0 can be calculated according to the formula 

TO = \ * J* 

A serious disadvantage of this conventional method is that the calculated propagation time x 0 is not accu- 
rate enough, which is demonstrated with reference to fig. 5, which illustrates the shape of a correlation function 
23 which is a combination of a correlation function 24 for the navigation signal 5* as received via the line of 
sight and two correlation functions 25, 26 for the navigation signal 5' as received via reflection. The propagation 

45 times Xi and x 2 of the reflection peaks 25, 26 are larger than the propagation time x 0 of the line of sight peak 
24, and their respective maxima A1 and A2 are shown as being smaller than the maximum AO of the line of 
sight peak 24. The widths W1 and W2 of the reflection peaks 25, 26 are substantially equal to the width WO 
of the line of sight peak 24. 

In the conventional method, the correlation function 23 is treated as if no multipath is present, i.e. the cor- 

so relation function 23 is processed based on the assumption that it is an undistorted correlation function. Then, 
the conventional method may, for instance, result in l(t x ) and l(t x +At x ) satisfying the stop criterion for the error 

signal e, and the propagation time x x is calculated as x x = t x + -At x , which is obviousiy not equal to the prop- 
agation time x 0 of the line of sight peak 24. 
55 . In fact, the time x x calculated by the conventional method is not even an accurate estimate of the position 
of the maximum of the peak 23, because the peak 23 is not symmetrical in the presence of multipath. The 
conventional method can only be improved in this respect by decreasing At, in order to sample the peak 23 at 
higher positions on the peak. This may, for instance, result in l(t y ) and l(ty+Aty) satisfying the stop criterion for 
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the error signal e, and the propagation time x y is calculated as x y = t y + ^Aty, which is a closer estimate of 
the position of the maximum of the peak 23 but, however, still does not equal t 0 . 

In fact, even if the conventional method were modified such that the highest point of the peak 23 were 

5 calculated accurately, it is inherently impossible for the conventional method to yield the correct value of t 0 in 
situations as illustrated in f ig.5 because, as can clearly be seen in f ig.5, the maximum of the combination peak 
23 is shifted with respect to the maximum of the line of sight peak 24. 

It is noted that in f ig.5 the peaks are shown in their ideal form, i.e. with sharp tips. However, as noted above, 
in practice the tips will be curved, which will cause a further increase in the deviation of the calculated value 

10 of the propagation time with respect to t 0 . 

The present invention is based on mathematical insight in the actual shape of the navigation signal 5' and 
the correlation function 23 in the presence of multipath. Therefore, according to one aspect of the invention 
the influence of multipath is eliminated, or at least minimized, by measuring the complete shape of the corre- 
lation function 23 of the received navigation signal 5', calculating an estimate of the multipath-induced contri- 

15 bution and the line of sight contribution to said correlation function, and calculating the propagation time x 0 
from the thus calculated line of sight contribution to said correlation function. 

This will be explained in the following discussion, wherein p(t) is a function representing the navigation 
signal 5, which corresponds to the reference signal as stored in the memory of the receiver 10. It is assumed 
that there are M ac reflection signals (multipath) contributing to the correlation function 23 of the received nav- 

20 igation signal 5', and that each of said M ac reflection signals can be represented by the function x,(t) = a . p(t- 
t,) - cos(ttt + 9,), wherein 
ai is the amplitude of a signal with index i; 
Ti is the delay time of the signal with index i; 
9| is the phase of the signal with index i; 

25 <q is the angular frequency of the carrier signal of the navigation signal. 

Then, in the presence of multipath propagation, the actually received navigation signal 5' as received by 
the receiver 10 can be written as 
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r(t> - £ xi(t) + n<t) « ]T ai • p (t-xj.) • cos (<Dt+9i> *n(t) (1) 
i-0 i=0 



Herein, the index i = 0 denotes the line of sight signal, while the index 1 = 1.2. ... M ac denotes the multipath 
35 signals. Further, n(t) denotes a noise component. 

In practice, the signal x,(t) is also multiplied by a data signal, which can be removed by multiplying the signal 
with an estimated data signal, as will be evident to a person skilled in the art since such data signal must also 
be removed in conventional coherent delay lock loop systems. For simplicity of notation, the data signal is left 
out of formula (1). 

40 It is observed that a,, x, and 9( are unknown variables which are representative for the signal with index i, 

and which may vary in time. According to maximum likelihood estimation theory, the best estimates a,, x, and 
9| of the unknown parameters a,, x, and 9, are those values which minimize L[s(t)] which is defined as 



L[s(t) ] - J [r(t) - 3(t) ] 2 dt (2) 



wherein s(t) is defined as 

Mac ^ac 

s(t) = X *i <t) = X ai • p(t-Ti ) • cos (cot + §i) (3) 

i-0 i=»0 

It is noted that in practice, where the signals are sampled, the above integral operation can be replaced by a 
summation over all samples in a certain time interval. 

A criterion equivalent to the above is that all partial derivatives of L(s(t)) are zero. It can be shown that this 
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criterion is fulfilled with 



max^ 

T 



f Mac . ^ 

Rej R x (t) - X a m - R(x-T ra ) • exp ( j§ m ) 

nr=0 



exp(- j9i) 



jo 



15 



20 



ai = Re 



f 



©i = arg 



Rx(ti) 



I 



a m - R(Ji-t m ) • exp( je m > 



Mac 

X a m - R(ti-x ra ) • exp( j§ m ) 
m=0 



\ 
) 



exp (- j8i) 



(4) 



wherein T{X} is defined as that value of x for which X(t) has its global maximum; 



25 



R x (t) 



p ( t-t) • exp (- jcot ) dt 
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is the in-phase (real part) and quadrature (imaginary part) downconverted correlation function, corresponding 
to the combination correlation function 23 discussed in the above example; and R(x) is the in-phase (real part) 
and quadrature (imaginary part) reference correlation function, corresponding to the correlation function 8 dis- 
cussed in the above example, defined as the correlation function as will occur in practice in the receiver, and 
normalized such that at the peak the maximum value is one and the phase is zero, such as illustrated in f ig.3F. 
According to the invention it is preferred that the reference correlation function is measured in the absence of 
noise and multipath, using a signal simulator, and that this shape is stored in a memory of the receiver as the 
reference correlation function, for instance in the form of a table or a function representing said shape. 

It will be clear that the objective of the invention is attained as soon as the equations (4) are solved, be- 
cause then the delay time t 0 of the line of sight signal can be considered calculated accurately by x^. However, 
the equations (4) represent recurrent relations between the optimal parameters %>, a 0 , 9 0 , x 1v d 1v 9 1f .... x Mac , 
&Mact 9^, for a specific value of M ac . Although in theory it will be possible to solve equation (3) directly by 
using an iterative matrix calculation, this approach costs so much calculating time that a real time implemen- 
tation is only possible with extremely fast and therefore large and expensive computers systems, if at all. 

In a more practical procedure according to the invention, the equations (4) are calculated iteratively, as 
described as follows with reference to f ig.6. 

In a first step 101, afirst peak correlation function Ro(t) is approximated as being equal to R x (x), according 
to the formula Ro(t) = R x (t). 

In a second step 102, parameters f 0 , a 0> 9 0 are calculated, in a manner to be described later, as optimum 
values for the largest peak in Ro(t). 

In a third step 103, the contribution of.said calculated parameters x 0 , a 0 , 9 0 is subtracted from R x {x) to yield 
a second peak correlation function R^x), according to the formula R^x) = R x (x) - a 0 • R(x-x 0 ) • exp(j9 0 ). 

In a fourth step 104, parameters x 1f § lf are calculated as optimum values for the largest peak in R 1 (x). 

In a fifth step 105, the contribution of said calculated parameters x 1f d 1v is subtracted from R x (x) to yield 
a new approximation of the first peak correlation function Ro(x) according to the formula Rq(x) = R x (x) - aV 

In a sixth step 106, new parameters x 0 , a 0 , 9 0 are calculated as optimum values for the largest peak in the 
new approximation Ro(x) as calculated in the fifth step. 

The third, fourth, fifth, and sixth steps are then repeated until a suitable stop criterion 201 is fulfilled. As 
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an example for such a suitable stop criterion 201, said steps may be repeated for a predetermined number of 
times, for instance 10. However, it is also possible to repeat said steps until the parameters x 0 , a 0 , 0 O and x 1f 
a 1f 0, have converged according to predetermined convergence criteria, as will be clear to a person skilled in 
the art. A convergence criterion which has proven to yield good results, is detecting whether the estimated delay 
x 0 and/or changes less than 0.1 ns between two successive iteration steps. 

At this point, for two signals the estimated contributions x 0 (t) and x,(t) to the received navigation signal 
5* are provided as Xo(t) = § 0 • p(t-x 0 ) • cos(©t + 9 0 ) and x\(t) = a^ • pCt-x,) • cos (a>t + 9,), and the estimated con- 
tributions Ro(t) and R,(x) tO A the in-phase/quadrature downconverted correlation function R x (x) are provided 
as Ro(t) = a 0 ■ R(x-x 0 ) • exp(j0 o ) and R^x) = 3, • Rfx-xO ■ expO^). 

In case an estimated contribution for a third signal is to be calculated, the contribution of said calculated 
parameters x 0 , a 0 , 8 0 and x 1f a 1f is, in a seventh step 107, subtracted from R x (t) to yield a third peak correlation 
function R 2 (t), according to the formula R 2 (x) a = R x (x) - a 0 ■ R(x-x 0 ) ■ exp(j0 o ) - a, - R(x-x0 • exp(j0i). 

in an eighth step 108, parameters x 2 , a 2 , 9 2 are calculated as optimum values for the largest peak in R 2 (x). 

In a ninth step 109, new parameters x 0 , § 0> 9 0 are calculated as optimum values for the largest peak in 
Roto = Rxto - St ■ R(x-X!) • expQeo - § 2 • R(x-X2) • expG8 2 ). 

In a tenth step 110, new parameters x\, a 1f 81 are A calculated as optimum values for the largest peak in 
Rito = R X M - a 0 ■ R(x-x 0 ) - exp(je 0 ) - a 2 • R(x-x 2 ) • expG'0^, using the parameters x 0 , a 0 , 9 0 as calculated in the 
ninth step. 

In an eleventh step 111, new parameters x 2 , 3 2 , 9 2 are A calculated as optimum values for the largest peak 
in R 2 (x) = R x (x) - a 0 - R(x-x 0 ) ■ exp(j9 0 ) - &y Rfc-^) • expg^), using the parameters x 0 , a 0 , % as calculated in 
the ninth step and the parameters x 1f S 1f 0, as calculated in the tenth step. 

The ninth, tenth, and eleventh steps are then repeated until a stop criterion 202 is fulfilled, which may be 
equivalent to the stop criterion 201. 

In an analogous way, the contribution of a fourth peak, and new values for the first three peaks, can be 
calculated iteratively; then, the contribution of a fifth peak, and new values for the first four peaks, can be cal- 
culated iteratively; etc., as will be clear to a skilled person. 

Stated differently, the contributions of all signals are not calculated simultaneously from the outset, but 
firstly the contribution of only one signal is calculated, then the contribution of a second signal is added and 
the contributions of both signals are optimized, then the contribution of a third signal is added and the contri- 
butions of all three signals are optimized, etc. It is possible to continue this procedure until the contributions 
of all M ac +1 signals are calculated However, it is noted that in practice the number M ac of multipath signals 
actually occurring is not known in advance and is even not necessarily constant. Therefore, after each cycle 
in the above-described procedure when the contributions of M+1 signals are calculated, wherein M is the con- 
sidered number of multipath contributions, the multipath error is not eliminated sufficiently if M is less than 
M ac , whereas an increased noise error is induced if M is larger than M ac . Therefore, in order to eliminate this 
problem, according to the invention it is preferred to solve the equations (4) successively for M = 0, 1, 2, .... 
Mmax. as explained above, and to stop this iteration process as soon as a suitable stop criterion is fulfilled, or 
when M=M max , whichever occurs first, wherein M max is a predetermined fixed value serving as an upper limit 
for the iteration process in order to prevent the calculation time from becoming too large. A suitable value is 
for instance = 3 or = 4. 

It is noted that, in the procedure as described in the above, the first peak to be calculated is referred to 
with index 0, whereas previously the index 0 was used to refer to the line of sight peak which, however, does 
not necessarily correspond to the firstly calculated peak. In the following, the index 0 will again be used to 
refer to the line of sight peak. It is noted that, after the above-described process is completed, the set of para- 
meters corresponding to the line of sight peak can be recognized easily because this will be the set having 
the smallest value of x,. 

Said iteration process may be stopped as soon as it is detected that the error increases, as expressed by 
the suitable stop criterion of formula (5): 

SRR(M) ^ SRR(M - 1) (5) 
Then, the parameter values obtained for M-1 are chosen as optimum values. 

In formula (5), SRR(M) is defined as the Signal-to- Residual-Ratio according to SRR(M) = a 0 2 /V(M) (6) 
wherein V(M) is defined as the variance of the residuals after estimation of ail parameters according to formula 
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fo wherein D is the delay window, i.e. the delay interval which contains the line of sight correlation peak and the 
peaks of the interfering multipath signals. In practice, a suitable length for D is twice the duration of a chip, 
and a suitable position for D is such that the estimated line of sight correlation peak is positioned substantially 
in the centre of D or within the first half of D, because multipath contributions are only expected at delay values 
larger than the delay value of the line of sight correlation peak. 
15 It is noted that in practice the in-phase/quadrature downconverted correlation function R x (t) is sampled, 

so that the integral operation of formula (7) can be replaced by a summation operation. 
As an alternative stop criterion, it is possible to use formula (8) instead of formula (5) 

SRR(M) ^ 5 (8) 

wherein 8 is a suitable treshold value which, for instance, may be chosen equal to that specific Signal-to-Noise- 
20 Ratio SNR which just leads to an acceptable error. 

Another alternative stop criterion is expressed by formula (9) 

wherein max( I res I ) is defined as the maximum absolute value of the residual function res in the interval 
25 D t res being defined as 



res(t) = R x (x) - £ a m - R(x-? m ) - exp ( j§ m ) (10) 

m-0 



In a simple implementation of the present invention, which has proved to yield sufficiently good results in 
most practical circumstances, M max is chosen as a fixed number equal to one. In that case, approximate para- 
meter values for the two strongest peaks, i.e. the line of sight peak and the strongest reflection peak, are cal- 
35 culated. Only in those cases wherein two or more reflection peaks of almost identical intensity are present will 
the remaining error be of practical significance, then. 

In the following, a description is given of the most preferred procedure for calculating parameters x, 3, 0 
as optimum values for the largest peak in a correlation function Rj(x), such as described with reference to, for 
instance, step 102. 

40 Firstly, it is determined which one of the complex samples of said correlation function provides the largest 

value for [Re(R|(x))] 2 + [lm(R((x))] 2 . The corresponding delay will be referred to as x max . In other words: 

tmax = max{ [Re(Ri(x) ) ] 2 - [ Im (Ri (t) ) ] 2 } (11) 

45 t 

A 

From said one sample, the phase 8 is estimated according to formula (12) : 

'-"-( gaagii) ««.«Kw»o ] 



6 jc/2 if Re (Ri (T max ) ) =0 and Im(Rt (tmax) ) >0 

§ = -Ti/2 if Re (Ri (x max ) ) -0 and Im ( Ri (T max > > <0 
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Once the phase is known, it is used to perform a phase rotation on the input signal in order to obtain a 
coherently demodulated signal according to formula ^13): 

R\{x) = Re(R,(T))cos0 + lm(R,(T))-sin9 (13) 
The thus obtained phase demodulated signal is used to estimate the delay and the amplitude. 
The delay t of a peak can be estimated in various ways. In the following, two basically different approaches 
will be discussed. 

Firstly, it is possible to use all samples of the peak, and to calculate interpolation values to the samples, 
in order to obtain a more or less continuous representation of the shape of the peak. Subsequently, the position 
of the maximum can be calculated using well known analytical methods like the Newton-Raphson method. This 
method provides an optimal delay estimate with respect to the noise variance, if the sampling interval satif ies 
the Nyquist criterion. 

Secondly, it is possible to use only a small number of samples of the peak. This method is based on the 
fact that the shape of the peak is known in advance, and is preferred according to the invention because it is 
much faster than the previously described method. Said number of samples can be as small as two, as will be 
15 discussed in the following with reference to fig.7. 

Fig.7A shows a correlation function, comparable to the correlation function 8 as shown in f ig.4. The delay 
of the peak is indicated as t 0 . Samples A and B of this function are shown, taken at two different delays t A and 
t 2 , spaced by the sampling interval At, such that t 2 = t, + At. The measured correlation magnitudes R'M) and 
R'i(t 2 ) are indicated as R1 and R2, respectively. 
20 In the following, the delay lying exactly in the middle of t, and t 2 will be indicated as t c , such that = t c 

- ^Ax and t 2 = t c + ^At. 

Further, the distance from t<. to t 0 will be indicated as t x , such that = t 0 - x x - ^Ax and t 2 = t 0 - x x 

25 + ~Ax. It will be evident that the measured correlation magnitudes R'|(x c - |ax) and R*i(x c + |ax) can be con- 
sidered as functions of t x , since t 0 and At are constant. 

With said measured correlation magnitudes R',(t c - ±Ax) and R\{t c + |ax), a shape parameter f is cal- 
30 culated, which is defined according to formula (14) : 

R(*o " t x + ^At) - R(t 0 - t x - ±Ax) 

f(x x ) = ^ \ (14) 

lR(*o - t x + ^Ax)| + |R(x 0 - x x - ±Ax)| 

35 Therefore, f is also a function of t x . In f ig.7B a graph is shown of f(T x ), the vertical axis having arbitrary 

units. It appears that f is a continuously descending function in a certain area around the maximum, which 
means that there exists a one-to-one relationship between each value of x x and each value off. The shape of 
the function f(x) depends on the shape of the correlation function, the characteristics of the receiver, and the 
magnitude of Ax. Since, as is mentioned above, the shape of the correlation function can be established in 

40 advance for each receiver, it is also possible to calculate the function f(x) in advance, to calculate the inverse 
function f- 1 , and to store this inverse function in a memory of the receiver, for instance as a table or a polynomial 
representation. 

During operation, two sample values R',(t c - |at) and R',(x c + ^Ax) of the peak are used to calculate f ac- 
45 cording to formula (14). Using the stored inverse function H, subsequently x x is calculated. 

A 

Then, t 0 is estimated as x 0 = t, + x x + -Ax = t c + x x . 

Preferably, t< and t^Ax are chosen as close to the peak as possible, more preferably at opposite sides of 
the peak, because this will yield the most accurate result. 
50 As an alternative, one might define f according to 

R(x c + Iax) 

f= f— (15) 

R(x c - ^Ax) 
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In another alternative embodiment, four sampled values are used, wherein f can be defined according to: 

R(*o -* x - f A*) - 3R(t 0 - x x - Iat) + 3R(t 0 - t x + ^At) - R(t 0 - t x + |at) 
f(Tj = ^ 2 2_ (16a) 
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or 

_ R(T 0 -t x - 2At) - 2R(t 0 -t x - At) 4- 2R(t 0 - t x + At) - R(x 0 - t x + 2At) 
H^xJ j^j (15D) 

Herein, N is a normalisation parameter, equal to the sum of several absolute measured values. For instance, 

5 N= I R(x 0 ) | orN = |R(t 0 -|at)| + |R(t 0 + ^At)| 

In theory, such a third order difference method offers the advantage of being less susceptible to multipath, 
because the differencing operations act like a decon volution, decreasing the pulse width and causing over- 
lapping pulses to become more separated. On the other hand, inaccuracies in the receiver, especially in the 
10 local reference signal, will affect the accuracy of the achieved results more seriously. 

As an alternative to formula (13), it is possible to use an incoherently demodulated signal according to, 
for instance, formula (13a) or (13b): 

R',(t) = [Re(R,(T))P + [lm(R,(T))]2 (13a) 



15 



20 



25 



35 



40 



45 



50 



R'Kt) = V[Re(R<T))P + Dm(R,(T))]2 (13b) 
The fastest way of estimating the amplitude a of the largest peak in a signal is by simply chos in g the largest 
sample value as being the estimate of the amplitude. The resulting error can be made acceptably small, for 
instance by chosing At arbitrarily small. It is also possible to use the delay t as calculated above to control the 
sampling means, such that sampling of the peak occurs exactly at (or very close to) the calculated top of the 
peak. 

To estimate the amplitude, it is also possible to use similar techniques as described in the above in relation 
to the delay. A normalized peak shape function g can be defined as 

gw = IS < 15) 

Using one measured amplitude sample x and an estimated delay t as calculated with the above method, 
the amplitude can be estimated according to: 

' a = ife (16) 

As an alternative, it is possible to use the second order derivative of the peak, according to: 
30 g(T j _ R(t 0 -t x -At) - 2R(t 0 -t x ) + R( To -t x + At) 

N 

and to proceed in a similar way as described in the above with respect to formula (15). 

The advantageous effect of the invention will now be illustrated by discussing an example whith reference 
to f ig.8. 

A standard GPS antenna was placed in an arbitrary position at the foot of a large building, and the GPS 
navigation signal was received during some time. This signal was processed in two receivers simultaneously, 
the first receiver being a standard GPS receiver operating according to the principle as discussed with refer- 
ence to f ig.4, whereas the other receiver was set up to operate in accordance with the principles of the present 
invention. was chosen to be equal to 1 (one). The iterations were stopped when x 0 changed less than 0.1 
ns between two successive iteration steps, or after 10 iteration steps, whichever occurred earlier. The number 
of (complex) samples was 20. The distance between the samples one-twentieth of the duration of a chip for 
the first ten samples, and one-tenth of the duration of a chip for the following ten samples. Formula (14) was 
used for f. The integration time was 1 s. 

Fig. 8 shows the fluctuations 8t of the delay time x 0 of the line of sight peak as calculated by both receivers, 
translated into a distance measure 8L according to the formula 8L = c8r, plotted against the time t. At time t=0, 
the calculated values were used as reference values and set at zero; the values plotted in f ig.8 for t>0 are the 
values as calculated and subtracted from said reference values. 

The solid line indicates the fluctuations occurring in the standard GPS receiver's calculations. Due to mul- 
tipath, the calculated distance measure fluctuates between extreme values which are about 60 m apart, with 
a standard deviation of the calculated distance measure of more than 10 m. 

The broken line indicates the fluctuations occurring in the second receiver set up to operate in accordance 
with the principles of the present invention. Although, probably due to noise, the calculated distance measure 
still fluctuates, the extreme values are only about 5 m apart, while the standard deviation of the errors is re- 
duced to 1.1 m, i.e. an improvement by a factor of 10. 

In the following, an exemplary embodiment of the processing device 30 according to the invention will be 
described in more detail with reference to fig. 9. 

As described above, a receiver 10 comprises an input (antenna) 11 for receiving the navigation signal 4. 
The signal 4 as received is led to an input 31 of the processing device 30, which input 31 is coupled to a first 
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input 35 of a multiplier 34. It is noted that, if desired, the signal 4 as received may be amplified and/or down- 
converted to a predetermined intermediate frequency, and then digitized, before being led to the first input 35 
of the multiplier 34, as will be clear to a person skilled in the art; for the sake of simplicity, the amplifying means, 
downconverting means and digitizing means are not shown in fig. 9. 

The processing device 30 comprises a generator 32 for generating a signal with frequency co. An output 
33 of the generator 32 is coupled to a second input 36 of the multiplier 34. In the multiplier 34, the signal 4 as 
received is demodulated, and the demodulated received navigation signal 5' is provided at an output 37 of the 
multiplier 34. 

It is noted that the signal to be processed is to be treated as a complex signal, for which reason the gen- 
erator 32 in fact has at least two outputs, one providing a signal which can be referred to as cos(ot), and the 
other providing a signal which has a phase difference of nJ2 with respect to the first signal and which can be 
referred to as sin(cot). The cos(©t) signal is used to process the real (in phase) part of the received signal, where- 
as the sin(o)t) signal is used to process the imaginary (quadrature) part of the received signal. Both processing 
operations are in principle identical. Therefore, for the sake of simplicity, in fig.9 only one output 33 of the gen- 
erator 32 Is Illustrated, and the combination of said output signals will in the following be referred to as exp(ja>t), 
as will be clear to a skilled person. 

The processing device 30 comprises a plurality of N multipliers 40 1f 40 2 , 40 3> 40 Nl referred to in the 

following as the multipliers 40 k , k being an integer in the interval from 1 to N, wherein N corresponds to the 
number of samples to be taken from the received navigation signal 5'. Each multiplier 40 k has a first input 41 k , 
a second input 42 k> and an output 43 k . Each first input 41 k of the multipliers 40 k is coupled to the output 37 of 
the multiplier 34 to receive the demodulated received navigation signal 5' r(t)-exp(j©t). 

The processing device 30 comprises a plurality of N integrator means 45 k . Each integrator means 45 k has 
an input 46k coupled to the output 43 k of an associated multiplier 40 kl and an output 47 k coupled to an input 
51 k of a digital signal processor 50, which may be a computer, as will be clear to a person skilled In the art. 
The integrator means 45 k are arranged for summation of the samples received at their inputs 46k during a pre- 
determined period T, delivering the result of this summation operation at their outputs 47 k , and then restarting 
a new summation operation from zero. 

The digital signal processor 50 has a clock input 52 coupled to an output 53 of an accurate clock means 
54. Further, the digital signal processor 50 has a control output 55 for controllably driving a local reference 
signal generator 60, which has an output 63 for providing the local reference signal 6. The local reference signal 
generator 60 has a reference input 61 coupled to a memory 62, in which information relating to the characteristic 
code of the navigation signal 5 is stored, as described above. Although this is not shown in fig.9, the memory 

62 may contain information relating to a plurality of characteristic codes of different satellites, while the local 
reference signal generator 60 may selectably receive information of only one of said codes, such as to select- 
ably tune to the navigation signal of a predetermined satellite, as will be clear to a person skilled in the art. 

The processing device 30 comprises a plurality of N-1 delay means 70 1t 70 2 , 70 3 , ..... 70m. 1f referred to in 
the following as the delay means 70 n , n being an integer in the interval from 1 to N-1. Each delay means 70 n 
has an input 71 n and an output 72 n , and is arranged to provide at its output 72 n a copy of a signal received at 
its input 71 n yet delayed over a delay time Dx n . In an embodiment, all delay means 70 n provide the same delay 
time D t . The delay means are coupled in series, such that the input 71 n of a delay means 70 n is coupled to the 
output 72 n . ! of a previous delay means 70„. The input 71 1 of the first delay means 70! is coupled to the output 

63 of the local reference signal generator 60, which is also coupled to the second input 42^ of the first multipl ier 
40^ The output 72 n of each delay means 70 n is coupled to the second Input 42^ of an associated multiplier 
40 n+1 . 

The digital signal processor 50 is programmed to drive the local reference signal generator 60 in accurate 
synchronization with the time periods as specified by the clock means 54. However, the digital signal processor 
50 may shift the beginning of the first chip C(1) of the Ideal reference signal 6 with respect to the beginning to 
of a time period over a certain bias shift t b , as will be explained later. Then, each multiplier 40 k receives at its 
second input 42 k a local reference signal 6 which is shifted with respect to the beginning to of a time period 
over a shift time x k according to the formula 

Jc-l 
n-1 

This shifted local reference signal 6 can be expressed as p k (t) = p(t-x k ). 

At its output 43 k , each multiplier 40 k provides a multiplied signal 7 which can be expressed as 
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r(t)exp(jcot)p(t-Tk). 

Therefore, each integrator means 45 k provides at its output 47 k a signal which can be expressed as 

T 

^ J r (t) • exp( j(ot) ■ p(t-tjc) 
t=o 

and which constitutes a sample at t = x k of the in-phase/quadrature downconverted correlation function R x (x). 

10 Said samples are "refreshed" after each integration period T, and read by the digital signal processor 50 

at its inputs 51 k . The digital signal processor 50 is programmed to calculate the sets of equations (4), preferably 
using the iteration steps as described with reference tofig.6 t In which said samples read at its inputs 51 k are 
utilized. Preferably, in the case of a GPS system, the Integration period T is chosen not to exceed 10 ms In 
order to get proper data estimates. However, this does not mean that it is absolutely necessary that the digital 

15 signal processor 50 calculates the entire sets of equations (4) once every 10 ms; where only ionospheric 
changes play a role of importance, it suffices to calculate the sets of equations (4) once every second or even 
less often. 

More particularly, in the case of a GPS system, it is possible to estimate the phase of the carrier wave 
once every 10 ms in a conventional way. From this estimated phase, the data and the frequency deviation of 

20 the. local generator can be determined. The frequency deviation thus determined is used for controlling the 
frequency of the signal generator and the code generator, such that their frequency error remains substantially 
zero. The calculated data, which is present in the form of phase jumps, is used to remove the data from the 
correlation samples by multiplying the correlation samples with the estimated data bits, after which the thus 
multiplied correlation samples can be averaged over 1 sec or more. In order to obtain precise estimates of the 

25 phase and the delay, it then suffices to calculate the sets of equations (4) once every second or even less 
often. 

Based on the calculated phase and delay estimates, the digital signal processor 50 controls the generator 
32 and the local reference signal generator 60, which comprise preferably a numerically controlled oscillator. 
More specifically, the generator 32 is controlled such that 0 O is approximately zero. By amending the bias shift 

30 t b , the local reference signal generator 60 is controlled in such a way that the peak is kept within the delay 
window D. Alternatively, it is possible to control the local reference signal generator 60 in such a way that the 
phase of the maximum of the correlation function of the received signal is zero. For more information, reference 
is made to Spilker, who describes a tracking loop and delay lock loop which are designed for assuring that one 
of the samples substantially coincides with the top of the estimated line of sight peak, while two other samples 

35 are situated symetricaliy around said peak. 

The digital signal processor 50 further has a data output 59 for feeding the calculated estimates to a suit- 
able display means such as a plotter, a printer or a CRT, and/or to a suitable storage means such as a computer 
memory, as is well known in the art. In practice, the results are fed into a computer, which also receives the 
estimated results for several other transmitters in order to be able to calculate the exact position of the receiver. 

40 This calculation is identical to the calculation performed in existing receivers. 

For completeness sake, it is noted that in reality two multipliers 34 are present, the first one receiving said 
signal cos(at) and the second one receiving said signal sin(cot) at their respective second inputs 36, such that 
their respective outputs 37 provide the in-phase navigation signal 5' and the quadrature navigation signal 5\ 
respectively. Further, it is noted that the circuitry of the means 40, 45 and 70 is present in two-fold, for proc- 

45 essing said in-phase navigation signal 5' and said quadrature navigation signal 5', respectively. 

It will be evident that the described embodiments of the invention may be amended or modified by a person 
skilled in the art without departing from the scope of the invention or the scope of the claims. For instance, in 
cases where it suffices to only process the real (in phase) part of the received signal, formula (4) simplifies 
into formula (4A): 

50 



55 



12 



EP 0 654 677 A1 



Ti = max J 



Re 



7 



R X (t) - X am* R(T-T m ) 



m=0 



exp(-j9 0 ) 



ai = Re 



exp(- j6 0 ) 



§0 35 arg (Rx<*0) ) 



(4A) 



Further, it is possible to use other suitable functions for the shape parameters f and g as described, as 
long as the used function for f is a continuously descending or ascending function. 



Claims 

1. In a signal transmission system containing a transmitter and a receiver, a method for processing a signal 
within said receiver, comprising the following steps: 

receiving a signal; 
demodulating the received signal; 

obtaining samples R x (x) of the demodulated received signal at a plurality of times t; 
calculating an estimate of a multipath-induced contribution to said demodulated received signal, subtract- 
ing said estimated multipath-induced contribution from said demodulated received signal to obtain an es- 
timated line of sight contribution to said demodulated received signal, and calculating a propagation time 
t 0 for said signal propagating between said transmitter and said receiver, said propagation time t 0 is cal- 
culated from the thus calculated line of sight contribution to said demodulated received signal. 

2. In a signal transmission system containing a transmitter and a receiver, a method for processing a signal 
within said receiver, comprising the following steps: 

receiving a signal; 
demodulating the received signal; 
providing a reference signal; 

performing a correlation operation on the received signal with said reference signal in order to obtain a 
correlation function; 

obtaining samples R x (t) of the correlation function at a plurality of times t; 
calculating an estimate of a multipath-induced contribution to said correlation function; 
subtracting said estimated multipath-induced contribution from said correlation function to obtain an es- 
timated line of sight contribution to said correlation function; and 

calculating a propagation time x 0 for said signal propagating between said transmitter and said receiver, 
said propagation time t 0 is calculated from the thus calculated line or sight contribution to said correlation 
function. 

3. Method according to claim 2, wherein said estimate calculating step further comprises the step of calcu- 
lating parameters § h x, and 0, as follows: 
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f Mmax ^ 

i i = maxJ Re R x {t) - £ a m - R (x-t m ) • exp ( j 9 m ) exp(-j§i> 
T m-0 

LV m*i J 



( M max 
a 3 - Re| R x (ii) - X ^m* RUi-t*) - exp< j§ m ) 



exp (- j9i) 



9i = arg R x (xj.) - £ '*nr R ( ? rt m ) ■ exp ( j 0 m ) 

m-0 

^ m*i J 



wherein 



R X < T > = r j r(t>- p (t-t) • exp<- jot) dt 



wherein R(x) is the reference correlation function; 

wherein a h X| and 9| are estimated values for a i? x } and 0i, respectively. 

Method according to claim 2, wherein said estimate calculating step further comprises the step of calcu- 
lating parameters aj, x, and 9| as follows: 



max^ 
i 



Re 



-V 



m=»0 



J 



exp (- j8o) 



V 



ai = Re 



Rx<ti) " 



Mmax ^ 
X a m - R(Ti-tm) 
m-0 



exp (-j9o) 



wherein 



9o = arg(R x (x 0 )) 



T 

R x (t) = | J r(t)- p(t-T) • exp(-jwt) dt 
1 0 

wherein R(x) is the reference correlation function; 

wherein aj, X| and are estimated values for a ( , xj and 9 i( respectively. 

Method according to claim 3, wherein a reference correlation function is measured in the absence of noise 
and multipath, using a signal simulator, and wherein the shape of the measured reference correlation f unc- 
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tion is stored in a memory of the receiver, for instance in the form of a table or a function representing 
said shape. 

6. Method according to claim 2, wherein said estimate calculating step further comprises performing the fol- 
lowing steps: 

a) calculating parameters x 0 , § 0> 9 0 as optimum values for the largest peak in R x (x); 

b) subtracting the contribution of said calculated parameters x 0 , a 0 . &o from R x (x) to yield R^x) accord- 
ing to the formula R^x) = R x (x) - a 0 .R(x-x\>)exp(j8 0 ); 

c) calculating parameters x 1f a 1t 9i as optimum values for the largest peak in R^x); 

d) subtracting the contribution of said calculated parameters x 1t fi 1a 9, from R x (x) to yield a new approx- 
imation of Ro(x) according to the formula 

Roto = Rx(t) - a, ■ R(x - *,) ■ expCGO; 

e) calculating new parameters x 0 , 3 0 , % as optimum values for the largest peak in the new approximation 
Ro(x) as calculated in step (d); and 

f) repeating the steps (b), (c), (d), (e) until a suitable stop criterion is fulfilled. 

7. Method according to claim 6, wherein the steps (b), (c), (d), (e) are repeated for a predetermined number 
of times. 

8. Method according to claim 6 X wherein the steps (b), (c), (d), (e) are repeated until at least one of the para- 
meters x 0 , §o, e 0 and t\, S 1( 9, have converged according to predetermined convergence criteria. 

9. Method according to claim 6, wherein, in order to calculate an estimated contribution for a third signal, 
the following steps are performed: 

g) subtracting the contribution of the calculated parameters x 0 , S 0> 9 0 and x 1f a 1f 9, from R x (x) to yield 
R 2 (t), according to the formula 

R2W = R x (t) - % ■ R(x - x 0 ) ■ exp(j9 0 ) - S 1 • R(x - x^ • expO^); 

h) calculating parameters x 2 , 3 2 , 02 as optimum values for the largest peak in R 2 (x); 

i) calculating new parameters x 0 , £ 0 > % as optimum^ values for the largest peak in 

Roto = Rx(x) - S, ■ R(x - x\) ■ exp(j9i) - a 2 . R(x - x 2 ) • expqS^; 
j) calculating new parameters x 1( fl 1( 9, as optimum^ values for the largest peak in 

Rito = Rx(t) a - S 0 • R(x - £ 0 ) ■ exp(j9 0 ) - 3 2 ■ R(x - x 2 ) ■ expflty, 
using the parameters x 0 , 3 0 > % as calculated in step (i); 

k) calculating new parameters x 2 , § 2 , 9 2 as optimum values for the largest peak in 

R 2 (x) = a R x (x) a - § 0 • R(t - x 0 ) • exp(j9 0 ) - ^ • R(x - • exp(j0,), 
using the parameters x 0 , a^,, 9 0 as calculated in step (i) and the parameters x 1f a 1f Gi as calculated in 
step (j); and 

I) repeating steps (i), (j), (k), until a stop criterion is fulfilled. 

10. Method according to claim 2, wherein a phase 9 of a peak is estimated from one sample R,(x max ) of said 
peak, taken preferably but not necessarily at or close to a top of the peak, according to: 



»fRe(RKw))<0 

8 = nil 

'f Re(Ri(x max ))=0 and lm(R,(x max ))>0 

9 = -n/2 

»f Re(R,(x max ))=0 and lm(R,(x max ))<0 
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11. Method according to claim 2, further comprising the step of calculating a function f(x x ) in advance, wherein 
the inverse function f- 1 is calculated and stored in a memory of the receiver. 

12. Method according to claim 11, wherein the function f(tx) is defined as 

R(t 0 - x x + ^Ax) - R(x 0 - t x - Iax) 

fdx) = f \ 

|R(xo - t x + ^Ax)| + |R(x 0 - t x - ^Ax)| 

13. Method according to claim 11 , wherein the function f is defined as 

R(x c + ±Ax) 
f = - 2 



R(t c - |at) 

15 14. Method according to claim 11 , wherein the function f is defined as 

R(xo - t x - |at) - 3R(x 0 - x, - ±Ax) + 3 R(x 0 - x x + ^Ax) - R(x 0 - t x + §Ax) 

f(Tx) = 2 ^~Ti ~ — 



wherein N is a normalisation parameter. 

15. Method according to claim 11, wherein the function f is defined as 

r/ x _ R(x 0 -x x -2Ax)-2R(xq-x x - Ax) + 2R(x 0 -x x + Ax) - R(x 0 - x x + 2Ax) 

f(Tx) = ^ 

wherein N is a normalisation parameter. 

16. Method according to claim 14, wherein N is equal to the sum of several absolute measured values, for 
instance, 

N = I R(t 0 ) I orN = |R(t 0 -|at)| + |R(x 0 + iAx)|. 

17. Method according to claim 2, wherein said samples R x (x) are obtained simultaneously. 

18. Device for processing a signal, comprising: 
a digital signal processor; 

a generator for generating a signal having frequency ©, which signal can be expressed as exp(-j©t), the 
generator being controlled by said digital signal processor; 

first multiplier having a first input coupled to receive a navigation signal and a second input coupled to 
an output of the generator; 

a plurality of N second multipliers, each of said second multipliers having a first input coupled to an output 
of the first multiplier; 

a plurality of N integrator means, each having an input coupled to an output of an associated second mul- 
tiplier, and an output coupled to an input of the digital signal processor; 

a local reference signal generator coupled to a control output of the digital signal processor and having 
at least one output for providing a local reference signal; 

a plurality of N-1 delay means, each of said delay means having an input and an output, each output being 
connected to a second input of an associated second multiplier; 

wherein the input of a first one of said plurality of delay means is coupled to the output of the local reference 
signal generator and to the second input of a first one of said plurality of second multipliers, and wherein 
the input of another delay means is coupled to the output of a previous delay means. 



55 



16 



EP 0 654 677 A1 

3 




Fig. 1 




EP 0 654 677 A1 




18 



EP 0 654 677 A1 



M(«)f 

M(ti+At) 




tl \2 



ti+Ay t2+At 



Fig. 4 



I(ty)_ 



Ktx) 




tx+ 1/2 At x t x +At x 



Fig. 5 



19 



EP0 654 



677 A1 



101 





APPROXIMATE R0 (r) 




^102 




CALCULATE INITIAL 
VALUES a OF 

£0. ao, do 






I ► 


, ^103 




CALCULATE R1 (r) ] 








CALCULATE OPTIMUM 
VALUES OF 
rl, ai, 61 








CALCULATE NEW 
APPROXIMATION 
OF R0 (r) 








CALCULATE OPTIMUM 
VALUES OF 
CO. ao, 6o 


201-O 


YES 


^REPEAT^ 



CALCULATE R2 (r) 



107 



108 



CALCULATE INITIAL 
VALUES A OF 

r2, &2, 92 



109 



CALCULATE 

A A & 

r0, ao, 8o 



110 



CALCULATE 

A A £ 

r1, ai, 01 



I 



112 



CALCULATE 

A A & 

r2, a2, 92 




F/g. 6 



20 



EP 0 654 677 A1 





21 



EP 0 654 677 A1 




22 



EP 0 654 677 A1 



CO 
CD 



v 



oc 

O uj « 



UJ 

o 



u 



7 



CO 
W 

CM 

in\| 



in 



CO 

-J CC 

<o 

Z h- 

o < 

u. Z 
Ul UJ 



a 

Ocv 
5 <o 
UJ 

S 



o o o 

Q 0) O 
DC 

a 



o 
in 



QC 

o 

< _ 

QC m 

UJ 



CN 

in 



DC 
O 

< CM 

DC m 
O * 

UJ 



i- cm/ 
CO to 



DC 
UJ 

5:o 
I- «* 
-J 
D 

s 



a. 

p co 
5 



co 




in 



DC 
O 

< co 
DC in 

UJ 



CM CO J 
CO <o y 



DC 
UJ 

Is 



CM 
CM "T- 



7 



o 
»- 

< z 

ac m 

UJ 



CO 
CO 



1^ 



£E 
UJ 



CO 



3 
2 



CM CO 
CM t- 



DELAY 
70 










UJ N» 
Q 



Kco 

CM 



DC 
UJ 



z 

_ o 

_l 

s 



z 



CM 
CM 





>- 




< 










z 









\z 

CM 



Z 

CM 



in 

CO 



23 



EP 0 654 677 A1 




BEST AVAILABLE COPY 



EP 0 654 677 A1 



European Patent 
Office 



EUROPEAN SEARCH REPORT 



Application N amber 

EP 94 30 8632 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Citation of document with indication, where appropriate, 

of i— 1 — 



US-A-4 669 091 (NOSSEN) 

* abstract; figure 3 * 

EP-A-Q 552 975 (NOVATEL COMMUNICATIONS) 

* abstract; figure 1 * 

US-A-3 753 123 (CARPENTER ET AL) 

* abstract; figures 1,4 * 



Relevant 
to claim 



1,2,18 
1,2.18 
1,2,18 



The present search report has been drawn Dp for all claims 



THE HAGUE 



Z March 1995 



CLASSIFICATION OF THE 
APPLICATION (InkCL6) 



G01S5/14 
G01S1/04 
H04B7/005 



TECHNICAL FIELDS 
SEARCHED (IotCL6) 



G01S 
H04B 



Haffner, R 



CATEGORY OF CITED DOCUMENTS 

X : particularly relevant if taken aJoae 

Y : particularly relevant If combined witb another 

document of the same category 
A : technological background 
O : non-written disclosure 
P : Intermediate document 



T : theory or principle underlying the invention 
E : earlier patent document, bat published on, or 

after the filing date 
D : document cited in the application 
L : document dted for other reasons 

A : member of the same patent family, conea^ndnif 



25 



